Diffractive photon production at high momentum transfer in ep collisions 
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The diffractive photon production at large momentum transfer and large energies is a probe of 
the parton dynamics of the diffractive exchange. In this paper we revisit the leading order (LO) 
BFKL equation approach for this process and estimate, for the first time, the differential and total 
cross sections considering the next-to-leading order (NLO) corrections for the BFKL characteristic 
function. We obtain a reasonable agreement with the DESY HERA data. 
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I. INTRODUCTION 



The description of exclusive diffractive processes has 
been proposed as a probe of the Quantum Chromody- 
namics (QCD) dynamics in the high energy limit (For a 
recent reviews see, e.g. Ref. It is expected that the 
study of these processes provide insight into the parton 
dynamics of the diffractive exchange when a hard scale 
is present. In particular, the diffractive vector meson 
and photon production at large momentum transfer is 
expected to probe the QCD Pomeron, which is described 
by the Balitsky, Fadin, Kuraev, and Lipatov (BFKL) 
equation I n the last years, the HI and ZEUS col- 

laborations at DESY-HERA have measured the exclusive 
production of p, <p, J/ty and 7 with hadron dissociation 
and large values of \t\, the square of the four-momentum 
transferred across the associated rapidity gap. The ex- 
perimental data for vector meson production are quite 
well described in terms of the BFKL formalism at lead- 
ing order 0-[ll[. In contrast, f or p hoton production, the 
analysis presented in the Ref. [12| indicate that the this 
approach describe the energy dependence of the cross 
section but it is unable to describe its \t\ dependence. It 
is important to emphasize that in [12j the free param- 
eters present in the BFKL formalism at leading order 
has been constrained using the experimental data for the 
total cross section and the differential cross section was 
predicted and compared with the data. This procedure 
is the opposed from that used in Refs. @, 0] in order to 
describe the vector meson data, where the experimental 
data for the differential cross section for a fixed energy 
were used to constrain the free parameters and the en- 
ergy dependence of the total cross section was predicted 
by the LO BFKL equation. The possibility that the use 
of a different procedure allow us to obtain a better de- 
scription of the experimental data is our first motiva- 



tion to revisit the description of the diffractive photon 
production in the BFKL formalism. The second one is 
associated to the fact that in the last years the next-to- 
lcading order (NLO) corrections to the BFKL kernel were 
determined [lj, [l4| and the origin of the instabilities in 
the perturbative series was understand and the problem 
was solved using methods based on the combination of 
collinear and small- a; resummations ^ allow us 

to include, for the first time, the NLO corrections for 
the BFKL characteristic function in the analysis of the 
diffractive photon production at large- \t\. This is a first 
step in the direction of a full NLO calculation, since we 
are not including the NLO corrections to the impact fac- 
tor associated to the photon-photon transition. It has 



23( and its analytical form 
Thus, we are aware that 



'Electronic address: 
t Electronic address: 



barros@ufpel.edu.br 



werner . sauter @ uf pel . edu.br 



been calculated in Refs. [lj 
was recently presented in [24 [. 
our phenomenological study is a educated guess to a com- 
plete NLO calculation. However, we hope that this gives 
some enlightenments about the underlying QCD dynam- 
ics. 

The paper is organized as follows. In the next Section 
we present a brief discussion about the diffractive pho- 
ton production at large-i in ep collisions. In Section Mil 
we present the formalism for the calculation of the cross 
section is presented and the next-to-leading order correc- 
tions to the BFKL kernel are briefly revised, as well as, 
the different schemes used to improve the convergence of 
the perturbative series. In Section lTVl we present our pre- 
dictions and compare with the available data. Finally, a 
summary of our main conclusions is presented in Section 

El 



II. DIFFRACTIVE PHOTON PRODUCTION 

Diffractive processes such as ep — > eXY, where X and 
Y are hadronic systems of the dissociated photon and 
proton, respectively, have been studied extensively at 
HERA in the last decades (For a recent review see, e.g., 
Ref. [l|). One of the cleanest of the diffractive processes 
is that of diffractive photon production (X — 7), with the 
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FIG. 1: Diffractive production of a photon with high-t ex- 
change (hard Pomeron) in a virtual photon-proton scattering, 
producing a hadronic (Y) system. 



final state photon having a large transverse momentum 
and being well separated in rapidity from the hadronic 
system Y. The cross section for the process of electron - 
proton scattering measured in HERA is directly related 
with photon-proton scattering by the relation [12j ] . 

do(ep -> e + 7 + Y) = da(-/*p -> 7 + Y) ^ 



where W is the center-of-mass energy of the 7*p system, 
r(IV) is the Weizsacker- Williams approximation of the 
flux of photons produced by the electron [25[ , 
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with IY 2 = ys and Q 2 , 
0.01 GeV 2 . 

Although parton fragmentation, production of a re- 
solved photon and the Bethe-Heitler process contribute 
for this process [26| . the exchange of gluons in the t 
channel dominates over the other contributions at large 
s, which implies that the leading contribution in high 
energy regime comes from the Pomeron exchange (See 
Fig. [T]). In this regime, the cross section for the process 
photon-proton can be fully calculated in the factor im- 
pact representation and is given by the convolution of the 
hard partonic cross-section with the parton distribution 
function in the proton (27l - [3lj |. The final expression is 



da(j*p -> 7 + Y) 
dtdxj 

^2(qj(xj, \t\)+qj(xj, \t\)) 



Ol 

^(7*<Z^79), (3) 



where G(xj, \t\) and qj(xj, \t\) are the proton gluon and 
quark distribution functions, respectively, and da/dt is 
the differential cross section for the subprocess j*q — > •yq. 
The basic idea is that at large - \t\, the gluonic ladder 
probes a parton inside the target and the struck parton 



initiates a jet which carries the fraction Xj of the longi- 
tudinal momentum of the incoming hadron, being given 



by Xj 



-t/(—t + My — m 2 ), where My is the mass 



of the products of the target dissociation and m is the 
mass of the target. The minimum value of 
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calculated considering the experimental cuts on My. In 
what follows we consider a fixed value x mul = 10~ 2 . Ex- 
perimentally, in diffractive photon production at large - 
\t\ a backscattering photon is produced with a small an- 
gle in the detector and the transverse momentum of the 
photon is counterbalanced by the momentum of the jet 
scattered in the frontal region. The photon and the jet 
are separated by a large rapidity gap among them. 

The diffractive photon production can be modelled in 
the rest frame of the proton by a sequence of three events: 
first, the fluctuation of the incoming photon into a quark- 
antiquark pair a long distance from the proton target. 
Second, the pair interacts with the parton q via the ex- 
change of a color singlet state (the gluon ladder). Finally, 
the final qq pair annihilates into a real photon. In the in- 
teraction of the ladder with the parton in the proton, it 
transfers momentum and dissociates the proton, whose 
fragments, aftermost, hadronize. This process has been 
studied in Refs. |27H31f and its description is closely re- 
lated to the diffractive production of vector mesons at 
large \t\ discussed in detail in Refs. [8l-[TTj]. The main 
advantage present in photon production is that the theo- 
retical calculations are simplified by the absence of a vec- 
tor meson wave function, with the only non-perturbative 
part being the parton distribution functions of the pro- 
ton. However, the cross section is suppressed relative to 
that of vector meson production by the electromagnetic 
coupling of the qq pair to the final state photon. In the 
next section we present a brief review of the description 
of this process in the BFKL formalism. 



III. BFKL FORMALISM 

The energy and momentum transfer dependences of 
the cross section for the diffractive photon production is 
directly related to the description of the differential cross 
section for the partonic process 7*9 —> jq, which can be 
expressed by [H, H3] , 
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where the scattering amplitudes A.(+ t ±) depend of the 
polarization states of the photons (represented by the 
plus/minus signs) through the impact factors for the 
transition 7* — > 7. As the Pomeron couples directly to 
the partons in the proton, the scattering amplitude fac- 
torizes in the impact factor associated to the transition 
of the virtual photon into the real photon, the gluonic 
ladder and the coupling of the Pomeron with the parton 
in the proton. At leading order the amplitudes are given 
by 
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where a is fine structure constant; a s is the strong cou- 
pling constant (comes from the impact factor and will be 
kept fixed); e q is the quark charges; sq is a energy scale, 
chosen to be c\ ■ Q 2 + C2 ■ \t\ (ci and C2 are free parameters) 
as the meson production case Q (see discussion below); 
l>j(v) is the BFKL characteristic function; Q 2 is the pho- 
ton virtuality and T(z) is the gamma function. These 
results have been derived using the Mueller- Tang [32| 
prescription for the parton - Pomeron coupling and only 
have considered the contribution associated to the low- 
est conformal spin (to = 0). In Ref. [33| these results 
have been rederived and it was shown that the only one 
non-zero contribution for highest conformal spin, to = 2, 
is small. 

The BFKL characteristic function u)(v) is in general 
expressed as follows 

u(v) = u s x{l) (7) 

where a s = (N c a s )/n, N c is the number of colors and 
the function and 7 = 1/2 + iv. The function \ is given 
at leading order by [2j-l5| 

X LO (7)=2V.(l)-V(7)-^(l-7) (8) 

where ip(z) is the digamma function. This characteristic 
function has been used in Refs. 12 9143 ll in order to esti- 



mate the diffractive photon production. However, several 
shortcomings are present in a leading order calculation. 
Firstly, the energy scale sq is arbitrary, which implies 
that the absolute value to the total cross section is there- 
fore not predictable. Secondly, a s is not running at LO 
BFKL. Finally, the power growth with energy violates s- 
channel unitarity at large rapidities. Consequently, new 
physical effects should modify the LO BFKL equation at 
very large s, making the resulting amplitude unitary. A 
theoretical possibility to modify this behaviour in a way 
consistent with the unitarity is the idea of parton sat- 
uration [34j], where non- linear effects associated to high 
parton density are taken into account. Another possi- 
ble solution, which is expected to diminishes the energy 
growth of the total cross section, is the calculation of 
higher order corrections to the BFKL equation. After an 
effort of several years, the next-to-leading order (NLO) 
corrections were obtained [l3l . [l4T |. with the x function 
being given by 



X(7) = X L °(7)+«,X NL0 (7), (9) 



where the correction term 

x nlo is given by [H [3 



X NL °(7) = C X LU (l) + 7 W'{l) + - 7)] - 7 W7) + 4>{l ~ 7)] 
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withC = (4- 7T 2 + 5/3 /N c ) /12, (3 = (1L/V C - 2iV/)/3 is number of flavours, ^ n \z) is the poli gamma function, 
the leading coefficient of the QCD (3 function, Nf is the £(n) is the Riemann zeta-function and 
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The main problem associated to these NLO contributions 
is that are so large that the problem appears perturba- 
tively unstable. In particular, they imply negative cor- 
rections to the leading order kernel, resulting in a com- 
plex functional structure (pole positions) that gives an 
oscillating cross section [35(. Moreover, there are also 
problems associated to the choice of energy scale, the 
renormalization scheme and related ambiguities. 

An alternative to cure the highly unstable perturba- 
tive expansion of the BFKL kernel was proposed in Ref. 
[l5j , who realized that the large NLO corrections emerge 
from the collinearly enhanced physical contributions. A 
method, the cj-expansion, was then developed to rcsum 



collinear effects at all orders in a systematic way. The 
resulting improved BFKL equation was consistent with 
renormalization group requirements through matching to 
the DGLAP limit and resummation of spurious poles. In 
this approach the kernel is positive in a much larger re- 
gion which includes the experimentally accessible one. 
This approach was revisited in Ref. [l8[ obtaining an ex- 
pression for the collinearly improved BFKL kernel which 
does not mix longitudinal and transverse degrees of free- 
dom and reproduces very closely the results from [l5| . 
The characteristic function proposed in Ref. [HI, de- 
noted All-poles hereafter, is given by 
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Another alternative to solve the spurious singularities 
present in the original NLO kernel was proposed in Ref. 
[3. Differently from Refs. QJ, QJ], where the calcu- 
lations were performed by employing the modified min- 
imal subtraction scheme (MS) to regulate the ultravio- 
let divergences with arbitrary scale setting, in Ref. [l6| 
they propose to solve the energy scale ambiguity using 
the Brodsky-Lepage-Mackenzie (BLM) optimal scale set- 
ting [36[ and the momentum space subtraction (MOM) 
scheme of renormalization. In this approach, the BFKL 
characteristic function is given by 



« V1 = X LU (7) 



a M OM{Q 2 )N c 



«mom(Q 2 ) 



(14) 



where omom is the coupling constant in the MOM 
scheme, 



«mom = ct s 1-1 ^-Tmom (15) 

L 7T J 

with T being a function of number of colors, number 
of flavours and of a gauge parameter (See Ref. [l6[ for 
details) and 
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Moreover, the function Q is the BLM optimal scale, 
which is given by 



Q 2 {v) =Q 2 exp 



, (17) 



with q = -2 jg 1 dx\n(x)/(x 2 - x + 1) « 2.3439. Finally, 
f is the NLO coefficient of the characteristic function, 
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section for the process 7*p — > 7 V is given by 



FIG. 2: (color online), v dependence of the characteris- 
tic function obtained considering NLO, NLO+All-poles and 
NLO+BLM approaches. For comparison, the behaviour of 
the LO BFKL characteristic function is also presented (a s = 
0.23). 



with 



x ^ cos(f ln(a;)) 
(l + x)y/x 



— -Li 2 i 
6 



x) 



(19) 



where Li2 (a;) is the Euler dilogarithm or Spence function. 

In Fig. [5] we present a comparison between the dis- 
tinct BFKL characteristic functions discussed above. For 
comparison we also present the LO BFKL function. We 
have that the All-poles and BLM approaches regularize 
the behaviour of the original NLO BFKL characteristic 
function at small- v and predict similar behaviours for the 
function u>. In particular, All-poles approach predicts 
that w(0) = 0.1 while BLM one predicts w(0) = 0.18. 
As we will show in the next section, this small difference 
have important implications in the energy dependence of 
the total cross section for diffractive photon production. 

In the next section, we will use Eqs. ©, (fT2"j) and (fT4")l 
as input in Eq. ([5]) in order to calculate the differential 
and total cross-section in the BFKL formalism. 



IV. RESULTS 

In this section we present our results for the differen- 
tial and total cross section and compare our predictions 
with the HI experimental data 12j . The differential cross 



da 
~dt 



d.r 



dtdxi 



(20) 



where we assume that x m i n = 10~ 2 . Moreover, the to- 
tal cross section is obtained by integration of the Eq. 
([2"0"]l over the |i| - range given by 4.0 GeV 2 < \t\ < 
36.0 GeV 2 , according with HI data jl2| . We will use 
the MSTW2008LO parametrization of the parton distri- 
bution functions [33], but we have checked that a differ- 
ent parton parametrization modifies slightly the predic- 
tions. In our calculations we have two free parameters: 
the strong coupling constant and the energy scale param- 
eter so- In what follows we assume that a s = 0.23, which 
is very close to the value used in Refs. [6T4TT1]. It is impor- 
tant to emphasize that in our calculations a s is fixed in 
the impact factors for all approaches studied, since they 
have been calculated at leading order, and running in the 
characteristic function when we are considering the All- 
poles and BLM approaches. In contrast with the vector 
meson production cross section, which is proportional to 
OQ, in the photon production case, it isproportional to 



a s . Following our previous studies 



is pi 



we have that 



*o = c\-Q +C2-|f|. Considering that the HI experimental 
data is for Q 2 = 10~ 6 GeV 2 we will assume for simplicity 
that Q 2 « 0, which implies sq = C2 • \t\. Consequently, 
in our calculations we have only one free parameter: the 
coefficient C2, which will be constrained by the HI data 
for the differential cross section as in Refs. 0, 0|- 

In Fig. [3] we present a comparison between our pre- 
dictions for the differential cross section considering the 
different approaches for the characteristic function and 
the HI data obtained for Q 2 = 10~ 6 GeV 2 and W = 219 
GeV. We adjust the value of so hi order to describe the 
larger number of experimental points. For the LO BFKL 
approach we obtain that C2 = 1.0 allow us to describe the 
data for |t| < 15 GeV 2 . Similarly, the value C2 = 0.2 al- 
low us to describe the data in the same kinematical region 
using the BLM approach. We have verified that if we ad- 
just the value of So in order to describe the data for \t\ — 
25 GeV 2 , these two approaches are not able to describe 
any other data point. In contrast, using the All-poles 
approach we only describe the data for \t\ >10 GeV 2 
(C2 = 6.0 x 10 -5 ). If we adjust C2 in order to describe 
the data for \t\ = 6 GeV 2 (c 2 = 1.0 x 10" 5 ) the approach 
does not describe the experimental data for larger values 
of \t\. Consequently, we can conclude that anyone of the 
approaches allow us to describe all experimental data. 
This conclusion also is valid if we consider the value of 
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FIG. 3: (color online). Differential cross section for the 
diffractive photon production at large-f considering differen- 
tial models for the characteristic function. Experimental data 
from [12]. 



a s as a free parameter. In this case the normalization is 
modified, which implies that a different of C2 also should 
be used in order to describe the data, but we have verified 
the shape of the ^-dependence is not strongly modified. 

Having fixed the free parameter in our calculations, we 
now are able to predict the energy dependence of the to- 
tal cross section considering the different approaches for 
the BFKL characteristic function. In Fig. [4] we present 
our predictions. We can see that the LO BFKL approach 
predicts, as expected, the steeper energy growth of the 
total cross section. In contrast, the All-poles approach 
predicts the milder energy dependence. In comparison 
with the experimental data, the better description is ob- 
tained by the LO BFKL approach. The BLM approach 
is able to describe the data at high energy, while the All- 
poles one (c 2 = 6.0 x 10~ 5 ) describes the data at low 
energy. For comparison we also present the result ob- 
tained using the All-poles approach for c 2 = 1.0 x 10~ 5 , 
which predicts a larger value for the normalization of the 
total cross section. 

Some comments are in order before we summarize our 
main conclusions. First, the current experimental data 
for the diffractive photon production are scarce and have 
been obtained in a limited energy range. Certainly, more 
data and for larger values of energy should allow us to 
obtain more definitive conclusions about the underlying 
QCD dynamics. The second, and more important, com- 
ment is related to the fact that in our calculations we 
are using leading order impact factors. In a full NLO 
calculation one should also consider the NLO corrections 
to the impact factors. In the last few years, the real and 
virtual corrections which contributes at NLO has been es- 
timated [19 23] and recently analytical expressions have 
been proposed (24{. Numerical results were presented 
in Refs. [23L [38| and indicate that the NLO corrections 
tend to decrease the value of the impact factors. Con- 



FIG. 4: (color online). Total cross section for the diffractive 
photon production at large-t considering differential models 
for the characteristic function. Experimental data from [12J. 



sequently, our estimates for the total cross section using 
the NLO-BFKL characteristic functions should be con- 
sider an upper bound. It is important to emphasize that 
the procedure used in our analysis was used in Refs. (39i — 
Elj ] in order to estimate the the double production of 
vector mesons (7*7* —> VV) and in Refs. [H, |43[ to 
calculate the 7*7* total cross section. Moreover, in Refs. 
[44H46I ]. the production of Mueller-Navelet jets was in- 
vestigated considering the NLO BFKL kernel and LO 
impact factors and, more recently, the description of the 
F2 structure function has been performed using a simi- 
lar approach [U H|| . All these studies demonstrate that 
the study of the diffractive photon production modifying 
only the BFKL kernel is justified in order to a first esti- 
mate of the magnitude of the NLO contributions for this 
process. 



V. SUMMARY 

The description of the high energy limit of the Quan- 
tum Chromodynamics (QCD) is an important open ques- 
tion in the Standard Model. During the last decades 
several approaches were developed in order to improve 
our understanding from a fundamental perspective. In 
particular, after a huge theoretical effort, now we have 
available the NLO corrections for the BFKL character- 
istic function, which allow us to improve the analysis of 
the processes which are expected to probe the underly- 
ing QCD dynamics. In this paper we have studied, for 
the first time, the diffractive photon production consid- 
ering the NLO characteristic functions proposed in Refs. 
[THHH. Moreover, the LO BFKL approach for the pro- 
cess was revisited. The free parameters in our calcula- 
tions has been constrained considering the HI data for 
the differential cross section at fixed energy and after 
we have predicted the energy dependence of the total 
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cross section. We obtain a reasonable description of the 
scarce experimental data. New data in a large energy 
range should allow us to obtain more definitive conclu- 
sions about the QCD dynamics. Moreover, the inclusion 
of the NLO correction in the impact factor is an impor- 
tant step which we expect to perform in a next future. 
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